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Abstract 

The RPA long range correlations are known to play a significant role in understanding the 
depletion of single particle-hole states observed in (e, e') and (e, e'p) measurements. Here the 
Random Phase Approximation (RPA) theory, implemented using the DIS force is considered for 
the specific purpose of building correlated ground states and related one-body density matrix 
elements. These may be implemented and tested in a fully microscopic optical model for NA 
scattering off doubly-closed-shell nuclei. A method is presented to correct for the correlations 
overcounting inherent to the RPA formalism. One-body density matrix elements in the uncorrelated 
(i.e. Hartree-Fock) and correlated (i.e. RPA) ground states are then challenged in proton scattering 
studies based on the Melbourne microscopic optical model to highlight the role played by the RPA 
correlations. Effects of such correlations which deplete the nuclear matter at small radial distance 
(r < 2 fm) and enhance its surface region, are getting more and more sizeable as the incident energy 
increases. Illustrations are given for proton scattering observables measured up to 201 MeV for the 
16q^ ^'^Ca, ^^Ca and ^"^Pb target nuclei. Handling the RPA correlations systematically improves 
the agreement between scattering predictions and data for energies higher than 150 MeV. 

PACS numbers: 21.10.Gv, 24.10.Hi, 25.40. Cm, 25.40.Dn 



* Present adress: School of Physics, University of Melbourne, Victoria 3010, Australia, 
t Present adress: Lawrence Livermore National Laboratory, California, USA. 



1 



I. INTRODUCTION 



Our understanding of the many facets of the nuclear structure properties has been and 
still is relying on the picture of independent particles moving in a mean potential. This 
picture stands at the foundation of the shell model which nowadays serves routinely as the 
basis of nuclear structure calculations, and is implicit to the self-consistent mean-field (i.e. 
Hartree-Fock) description of nuclear ground states. For independent particle motion, the 
occupancy associated to nucleon orbitals is 1 or depending upon whether the single-particle 
level is below or above the Fermi energy, respectively. 

It is only recently that the quenching of shell model occupation probabilities has been 
disclosed in a dedicated series of experiments in which incident electrons serve to map 
detailed structure properties hard to reach using other probes. First hints revealing such a 
quenching came in measurements of electrons scattering from 206p]^ g^^^j ^°^T1, from which 
the 3s 1/2 proton radial wave function was determined. Its shape is peaked in the central 
region, and close to expectations for a 3si/2 wave function. Minor adjustment of the 3s-hole 
strength provided an improved data prediction jj]. Evidence for partial occupancy for this 
orbital was provided later on from a joint analysis of (e, e') and feje'p) experiments. The 
3si/2 orbital was found to be depleted by a (18 ± 9)% amount j2, |3|. Today the absolute 
occupation probability of this proton orbital is evaluated to be 0.76 ± 0.07 |^. 

Further detailed information on the single-particle structure have been recently gained 
through measurements of the spectral function S{E,k), where E and k are the removal 
energy and momentum, respectively, of a proton in (e, e'p) knockout experiments. For ^°^Pb, 
these measurements performed at high binding energy and momentum transfer show that 
mean-fieldpredictions are lying far below the data, highlighting the need for consideration 
of tensorJH 01 as well as short- 0, H, 0| and long-range correlations beyond the mean field 



|10l IllL I12L Il3j . A wealth of methods and models have been adopted to tackle this issue. 
These are the Green's functions metho d [l^ [l5| , the variational Monte-Carlo method 0, 



17l |. the correlated basis function theory [l8|, the particle- vibration model [19|, the dispersive 
optical model extrapolated to bound state region 0, H H 113, and the Random Phase 



Approximation 2ihe 2^ l27|, |28[. Among the correlations which have been considered 
so far, the long-range ones appear important for curing the deficiencies associated with the 
mean-field predictions. 

In the present work we investigate the impact that nuclear long-range correlations have 
on the interaction of nucleons incident on doubly-closed-shell nuclei, among which ^°^Pb, a 
nucleus for which many scattering observables have been measured. In the past, detailed ex- 
perimental information on nuclear structure as gained from electron scattering measurements 
played a key role in building effective NN forces and mass operators for nucleon scattering 
studies in the folding model framework (iol . l30| . Now, that a successful NA microscopic 
optical model (OM) based on a g- matrix interaction has been established in r-space |3l|, it 
is timely to push the limits of its predictive power using various microscopic picture infor- 
mation. Several studies along this line have already been published. For example, no core 
shell model wave functions have been adopted in successful interpretations of proton scat- 
tering measurements for ^^C and light nuclei below and at the neutron drip-line 3^, 3^ 3^ . 
Hartree-Fock predictions based on Skyrme forces have also been challenged in proton and 
neutron elastic scattering studies at medium energy to provide estimates of neutron skin 
thickness in ^°^Pb. Here, the correlated ground states of stable doubly-closed-shell nuclei, 
built using the finite range, density dependent DIS force (3^ in the self-consistent RPA 
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theory [36|, are used instead and thoroughly tested. 

Our paper is organized as follow. The main features of the fully antisymmetric, micro- 
scopic NA optical model are described in Secjlll Section ITTTl includes a brief presentation of 
the HF+RPA theory for establishing our notations, and describes the method used to fix 
the well-known double counting problem. RPA predictions are compared to experimental 
data for charge and neutron radial shapes of ^ospi^ its ground state. One body density 
matrix elements in the correlated ground state are then provided. Finally, optical model 
predictions based on HP and HP+RPA one body density matrix elements are compared in 
Sec JIVI to various scattering observables in the 40-201 MeV incident proton energy range for 
i6q^ ^°Ca, ^'^Ca and ^°^Pb, and to scattering predictions based on the Skyrme SkM* js^l 
force. 



II. MICROSCOPIC OPTICAL POTENTIAL FROM THE MELBOURNE g MA- 
TRIX 



The full details of the Melbourne g matrix optical potential may be found in Ref. j31[, 
to which we refer the reader. We present a brief summary of the derivation of the potential, 
highlighting those points relevant to the use of RPA densities in its calculation and the 
observables obtained therefrom. 

In folding models of the optical potential, one starts with a credible effective A^A^ inter- 
action. In the case of the Melbourne potential, the effective A^A^ interaction is the g matrix 
derived from the Bonn-B A^A^ interaction [s^ . The g matrix for infinite matter is a solution 
of the Bruckner-Bethe-Goldstone equation in momentum space, viz. 

g (q', q; K) = V (q', q) + / ^ (q', k') ^E^^^^^'^E%,K)f ^' ' 

where Q(k' , K; kf) is a Pauli operator and medium effects are included in the energy denom- 
inator. Effective g matrices are obtained in coordinate space for finite nuclei whose Fourier 
transforms best map those momentum space solutions. Those g matrices so obtained contain 
central, tensor, and two-body spin-orbit terms. They also are constructed over all two-body 
spin and isospin channels, allowing for a self-consistent specification of proton and neutron 
scattering, as well as charge exchange reactions. Those g matrices are then folded with the 
ground state density matrix elements to give the optical potential for elastic scattering. 
The optical potential derived therefrom can be cast in the form 

U (r, r'; E) = 6 {r - r') ^ / (^;(s)^z)(r, s; E)^p{s) ds 

a/3 

= t/z5(r; E)5 (r - r') + Ue{v, r'; E) , (2) 

where the subscripts D and E designate the direct and exchange contributions, respectively. 
Nuclear structure information enters in via the one-body matrix elements pa/i and in the 
specification of the bound state single particle wave functions and ipp. In terms of the 
RPA (or HP) ground state |0), that density matrix element is paj3 = (0|a^a^|0) (see SecJIII|). 
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The main source of non-locality in the optical potential is from the exchange term. The 
direct term resembles a gp-tjpe optical potential and by definition is local. The form of 
the exchange term necessarily does not follow this construction: the exchange terms in the 
folding require that the sum is over explicit effective NN two-body amplitudes. As such, 
direct comparisons are not possible between this form of the optical potential and those which 
are local, as constructed from nonlocal NN amplitudes through local approximations, or as 
specified phenomenologically as sums of Woods-Saxon form factors. 

To obtain the observables for scattering, the optical potential so obtained is used in the 
nonlocal integro-differential Schrodinger equation, viz. 



^(r) = j f/(r,r')^(r')t^r' , (3) 



where Vc{r) is the Coulomb potential, and the terms due to the intrinsic spin of the system 
have been suppressed for simplicity. The code DWBA98 [3^ is used to calculate the folding 
potential from the effective A^A^ g matrices and obtain the relevant scattering observables. 

At low energy, the averaging over the coupling to the nonelastic channels represented by 
the g matrix is no longer valid and the derivation of the optical potential must be done 
in terms of explicit channel coupling to open and closed channels. Such has recently been 
constructed in terms of the collective model |i4Q . i4ll | . 



III. NUCLEAR STRUCTURE 

As an introduction to this section, it is important to mention that our approach is not fully 
consistent. On the one hand, we use the g-matrix as an interaction between the projectile 
and the nucleons in the target, whereas on the other hand, to calculate nuclear structure, 
we consider effective interactions which have been separately adjusted. As long as we focus 
on studying medium energy scattering, one can find justifications for proceeding in this way. 
However, at low energy, this approach would be more questionable and it is likely that the 
derivation of an optical potential in the frame of a more fundamental theory (as in |4^ ) 
should be considered. 



A. The mean field approximation 

The simplest description of the nuclear structure is provided by the self consistent mean 
field theory, which is also called Hartree-Fock (HF). There, the ground state is a Slater 
determinant constructed with individual particle states which are solutions of the HF equa- 
tions. In this work , w e use the HF results obtained using two different interactions. One is 



the Skyrme SkM* 37[ interaction , and the other one is the finite range, density dependent 
DIS interaction [35]. The details of the HF formalism used with the DIS density dependent 
interaction can be found in 4^ 3 • 



In order to calculate the one-body matrix elements of SeclHl it is convenient to 
express the HF ground state in second quantization as 
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The above product contains only occupied states labeled "h" (hole states) according to the 
usual terminology. The creation operator associated with the creation of a hole in a HF 
single particle state is defined with: (ph{r) = {r\a'^\OHF)- 

By introducing these notations, the matrix elements p^^ read 

p^fS = {HF\a+a^\HF) . (5) 
and are diagonal {ph,h = 1 , Pp,p = 0) in the HF approximation. 



B. Description of the ground state beyond the HF approximation 

The density dependent effective interaction DIS has successfully been used in various 
extensions of the mean field theory. Among them, the one of interest for our study is 
the microscopic description of collective excitations for closed shell nuclei as described in 
reference js^. We recall some essential features of this approach and make the link with the 
usual RPA theory. This will permit us to define the two variants of correlated ground states 
that we propose for the description of the target. 



1. Ground state correlations induced by collective excitations 



The approach of 3^| is based on the quadratic form introduced a long time ago to 
study the stability conditions of the HF solutions. It is obtained by performing a Taylor 
expansion of the energy E up to second order in the variation of the density matrix around 
the equilibrium HF density (p^^^)- The quadratic form in question is expressed in terms of 
the matrix 



A B 

B* A* 



with elements 



and 



^{ph),(p'h') — ^pp'^h,h'{^p 



B. 



{ph),{p'h'} 



dpp'h' 



p=p 



(0) 



dpphdp. 



h'p' 



(6) 
(7) 
(8) 



p=p 



(0) 



where and eh are the HF single particle energies for a particle state and a hole state, 
respectively. This matrix is used to define a set of RPA equations 3^ 4^, namely 



A B 

B* A* 



X 
Y 



LU 



X 
-Y 



(9) 



where u is a. set of eigen-values corresponding to a set of eigen-vectors with components 
X and Y. The definition of the matrix (jH)) presents the advantage to show that, due 
to its explicit dependence on the density, the particle-hole matrix elements of DIS must 
contain the so-called rearrangement terms besides the usual ones. Notice also, that one 
retrieves the usual particle-hole matrix elements when the interaction does not depend on 
the density. Once such prescription is adopted for defining the particle-hole vertices, the 
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ap pro ach developed in j36| follows closely the standard RPA theory as described extensively 
in j45| . Below we only give the relevant definitions that introduce the quantities of interest 
for this work. We express the formalism in a representation that accounts for rotational 
invariance and refiection symmetries of the nuclear interaction and the mean-field as well 
(see Appendix^. Creation and annihilation operators are defined through a Bogolyubov 
transformation 

p,h 

e,(.,.,M) = KmKa)^^^ M) + X-;_,)A(,,,)(7r, J, M) , 

p,h 

which mixes the creation and destruction operators, A^^^(7r, J, M) and y4(p_/i)(7r, J, M) re- 
spectively, of independent particle-hole pairs with definite angular momentum and parity. 
The amplitudes X and Y are the components of the solutions of the RPA equations defined 
in (jH)). Since we work within the quasi-boson approximation, the Bogolyubov transformation 
is nothing but a canonical transformation between two sets of bosons. Excitation modes 
of the nucleus are then defined through the action of any creation operator 0+ onto the 
quasi-boson vacuum \0) of the destruction operator O. This is expressed as follows 

K,(vr,J,M)) = e+^^,^,,)|0) , (11) 

0«,(.,j,M)|O) = V2,7r, J,M . 

The quasi-boson vacuum can be constructed explicitly from the vacuum \HF) of the 
^ip,h) ^) operators. According to j45| it reads 



J 



|0) = Ne^\HF) , (12) 

with 

Z = 

2 

7r,J {ph),(p'h') 

and the normalization defined as 

= {HF\0) . 

This form shows clearly that the quasi-boson vacuum is a superposition of (2p-2h), 2 (2p-2h) 
... n (2p-2h) excitations coupled to zero angular momentum as it should, since the total 
spin of the ground state is zero for the nuclei under consideration. In the present work and 
for future applications to inelastic scattering we assume that the quasi-boson vacuum (fT^ 
and the excited modes (fTT|) provide a reasonable description of the ground state and nuclear 
excitations of the target. 

At this stage it is worth pointing out that there exists another explicit form of the 
correlated ground state that has been derived 12^] by summing up the RPA diagram to all 
orders. This important work shows that the resulting ground state, denoted here as \RPA), 
has exactly the same structure as the quasi-boson vacuum, but it reveals also that the quasi- 
boson counts twice the lowest order term of the perturbation theory. How it affects mean 
values of one body operator is now shown on the matrix elements of the one-body density 
operator. 
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2. One-body density matrix for the RPA ground state 



The one-body density matrix calculated in correlated ground states is no longer diagonal 
but contains all the elements of the form ph ^' and Pp^p'. The non-diagonal particle- hole 
matrix elements vanish because of the structure of the ground state. Besides, on the account 
of symmetries it can be shown that the density matrix reduces to diagonal block matrices 
labeled by r) and independent of the projection m of the angular momentum j, namely 

P{a),{P) ~ ^1^,1/3 ^jajp ^TcT/s P(na,la,ja,Ta),iniS,la,ja,ra) ■ (l^) 

Finally, it is often convenient in the formalism to perform the summation over m in advance 
and to consider the following quantities instead 

P(a),(/3) = ^P{a),{P) = (2ja + l)P(a),{/3) • (14) 
m 

We next provide expressions for these quantities in the cases of the quasi-boson vacuum 
and RPA vacuum 

m m 

The calculation in the quasi-boson vacuum is straightforward. We only give the result 
for the particle and hole cases, respectively, as 

i,J,iT,h 

and 



i,J,iT,h 



with the definition 5(a), (/?) = ^i^.i^^ic j^'^r^r^- 

In order to calculate the RPA one-body matrix elements one refers oneself to |2J| where 
expressions of the occupation probabilities of single particle orbital in the RPA state can be 
found. Although such probabilities involve only diagonal matrix elements of the density, it 
is not difficult to generalize an expression for the non-diagonal ones. It turns out that the 
one-body matrix elements in the RPA state and those in the quasi-boson vacuum ()16p differ 
only by the lowest order contribution in the perturbation theory. The correction terms are 
given for particle and hole cases, respectively, as 

Ap(a),(/3) = --()(a),(/3) 2^(2J + 1) 2^ 7- — — T , (17) 

pth^,h ^^(P'^h') + ^(c',h))y^(p',h') + ^W,h)) 

and 



j,^ p^pfh, i^ip',h') + e{p,«)) (e(p',fe') + e(p,/3)) 
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where the e^p^h) = — are the free particle-hole pair energies, and ^,-) the values 

defined in (jH} for particle-hole pairs with good angular momentum J and parity vr. With 
these notations, the RPA density matrix reads 

phm = P("),(/3) + ^p{c)m ■ (18) 

This expression is folded with the Melbourne g-matrix (see Q) and the optical potential so 
obtained is then used to calculate elastic scattering observables. 



3. Structure of correlated ground states 

From inspection of the vacuum structure (jl2j) as outlined in Appendix ^ it is clear that 
the 6a amplitudes (see IB3|) provide a direct measure of ground state correlations. Taking 
into account the (2 J + l)-fold degeneracy of the 9a 's in each (tt, J) subspace, the ratio 

is a measure of the relative importance of each subspace, with Jfief taken as the multipolarity 
of the one which provides the main contribution to the overall correlations (here, J^ef = 
3). These ratios shown in FigQfor ^osp]^ indicate that some natural and unnatural parity 
states of all (vr, J) subspaces, even high spin ones, are worthy of consideration for building 
the correlated ground state. 

As the correlations are smearing out the occupation probability distribution of proton 
and neutron single-particle levels around their respective Fermi energies, the radial g.s. 
densities get depleted towards the nuclear center. This effect can be seen in FigO where 
measured charge and neutron distributions are shown together with our HF and HF+RPA 
predictions for ^°^Pb. Calculated root mean square (rms) radii of proton, charge and neutron 
distributions as well as neutrons skins are gathered in Table |l] for ^°^Pb as well as for ^^O, 
^°Ca and ^^Ca. A good overall agreement between the RPA predictions and experimental 
values is obtained. 



IV. ANALYSES OF SCATTERING OBSERVABLES 

In order to test the predictions of the OMP described above, an incident proton ex- 
perimental database was built, comprising differential cross sections cr{9)/aRuth, analyzing 
powers Ay (9) and spin rotation functions R{9) and Q{9). References to these data are pro- 
vided in Table ini only for ^°^Pb. The incident energies of present interest are limited to the 
40-201 MeV range where the Melbourne OMP is most successful [s^l- For all the compar- 
isons between model predictions and experimental data shown below the continuous and 
dashed curves represent the OMP calculations based on one-body density matrix elements 
of correlated (RPA) and uncorrelated ground states (HF), respectively. 



A. Incident protons 

Proton scattering experiments have provided a wealth of valuable information on angular 
distributions for various observables at many incident energies. For this reason, the proton 
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FIG. 1: Values of the quantities Oa'"^ defined in the text. We present the contributions for ^'^^Pb, 
a = 1 — > 20 first states of each (7r,J) block. 
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FIG. 2: Charge and neutron radial densities of ^^^Pb. Comparisons between experimental data 
4^ . 47 1 (dotted curves), correlated (full curves) and uncorrelated (dashed curves) calculations. 
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HF 


5.432 


5.475 


5.567 


+0.135 




HF+RPA 


5.467 


5.504 


5.592 


+0.125 



TABLE I: Proton, charge, and neutron rms radii for ^^O, ^'^Ca, ^^Ca and ^''^Pb. Comparisons 
between present HF and HF+RPA predictions, and experimental values . The neutron skin Ar^p 

are 



is defined as Ar„p = {r'^y/^ — {rpY/"^ . The estimated {r'^Y^'^ and Ar„p values of Refs. 4^, 
from systematics. 



Energy (MeV) Ref. Energy (MeV) Ref. 
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[52] 


104.4,121.2 


[53] 


45,47.3 
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156 




61.4 


m 


160 




65 


[58] 


182.4 


[53] 


79.9 


[53] 


185 


[59] 


97 


m 


201 





TABLE H: a{e)/ajiuth, Ay{0), R{e) and Q{e) database for proton scattering off ^ospb. 

database we have formed in Table |n] serves as the main play ground for detailed OMP 
analyses. We also show some illustrations for the three other stable doubly-closed-shell 
nuclei ^^O, ^"Ca and ^^Ca. 
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1. ^"^^Pb 



The differential cross sections to be discussed below are normalized to Rutherford scat- 
tering cross sections to magnify differences existing between our OMP predictions and scat- 
tering data. This comparison is shown in the upper panel of FiglHl Similar comparisons for 
Ay [6) are shown in the lower panel of Fig|21 

For the comparison between solid (RPA-based) and dashed (HF-based) curves for cross 
sections, it turns out that the former is systematically lower over most scattering angles. 
Compared to OMP predictions based on the Hartree-Fock ground state density matrix, those 
using the RPA one are all in closer agreement with the spread of cross section data except 
maybe at lower incident energies where the calculated minima seem too deep. This is a 
known low energy shortcoming of the g-folding model which has been discussed previously 
(g^I- Nevertheless, the agreement between RPA-based calculations and measured differen- 
tial cross-sections is spectacular, especially at the higher incident energies, where HF- and 
RPA-based calculations differ the most. Extending the comparison from experimental cross 
sections to analyzing powers, it may be seen in the lower part of Fig|21 that the correlated 
ground state specifications lead to a excellent overall OMP description of the Ay[6) data 
spread, especially at medium angles for energies E > 150 MeV. 

A similar statement is made for the spin-rotation functions R{0) and Q{6) measured 
at 65 and 201 MeV, respectively. As can be seen in Figl^ the phasing and amplitude of 
these measured observables are well accounted for by our OMP calculations, although these 
observable predictions do not seem very sensitive to RPA correlations. 



2. Other doubly-magic nuclei 

Calculations were also performed for protons incident on the other stable doubly-magic 
nuclei ^^O, ^°Ca and ^^Ca. Although these calculations were performed for all incident 
energies where experimental data are available, FigElonly displays comparisons at highest 
energies, where the difference between HF- and HF+RPA- based OMPs is the most striking. 
Those results are representative of the agreement obtained over the 60-201 MeV range. For 
these doubly magic nuclei, comparison between calculations using correlated and uncorre- 
lated ground state density matrices, and the experimental data, allows us to confirm the 
conclusions of the p-|-^°^Pb scattering study made above with a larger data sample. 



B. Incident neutrons 



Although some neutron scattering data is available (HsI . 661 67l | at neutrons energies higher 



than 40 MeV, those data sets (with the notable exceptions of [66|, |67|) do not extend far 
enough in angles to allow for discrimination between the nuclear structure models used as 
a basis for our OMP analyses. Thus, those datasets can be described in a satisfactory way 
by our OMP using either HF or RPA one-body density matrix. Moreover, when comparing 
incident proton and incident neutron calculations, no effect specific to incident neutrons was 
observed, and like for incident protons, the RPA-based neutron-nucleus OMP calculations 
predict cross sections that are systematically lower at large angles than their HF counter- 
parts. Nevertheless, the scarcity of high energy, large angular range neutron scattering data, 
calls for new measurements of the quality of those in 6g, 6^, maybe at higher energy. 
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FIG. 3: Differential cross sections (T{9)/ajiuth and analyzing powers Ay{9) for protons incident 
on ^'^^Pb. Comparison between data (symbols) and OMP predictions based on correlated (solid 
curves) and uncorrelated (dashed curves) descriptions of ground state. Cross sections are offset by 
factors 10, while analyzing powers are shifted by 2. 
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FIG. 4: Spin rotation functions K{0) and Q,{9) at 65 and 201 MeV for protons incident on ^'^^Pb. 
Comparisons between experimental data (symbols) and present OMP calculations for correlated 
(solid curves) and uncorrelated (dashed curves) ground state descriptions. 



We conclude this analyses with making the statement that the RPA correlations have 
sizeable impacts on the OMP predictions only at the higher incident energies of present 
interest and for center-of-mass scattering angles larger than typically 6 ~30°. This statement 
is relevant to ^^O, ''^Ca, ''^Ca and ^osp]^ target nuclei. 



C. Discussion 



1. Probing ground state correlations 

In Sec. IIV A| we have shown that similarly to electron scattering, nucleon scattering is 
sensitive to small details of the nuclear structure of the target nuclei, such as those stemming 
from the presence of long range correlations in the target ground state. Moreover, including 
such correlations do improve the agreement between calculated and measured scattering 
cross sections. Next comes the difficult question of identifying the features of the correlated 
density matrix that nucleon scattering is sensitive to. Looking at FigOl can provide us 
with hints to that effect: the differences between HF- and RPA-based calculations can be 
seen to be stronger at large angles, suggesting that such differences appear when more 
interior regions of the target are probed. Re-plotting the p+^°^Pb scattering cross sections 
as functions of the momentum transfer q (see FigE)) confirms that indeed, for all energies, 
differences between HF- and RPA-based calculations are associated with values of q larger 
than 1.7 fm~^, and thus deeper regions of the target. Fig |21 displays the radial charge density 
of 208p|-, calculated with (solid curve) and without (dashed curve) RPA correlations in the 
ground state, showing the well known effect of RPA correlations, i.e. depleting the interior 
of the density distributions and enlarging the distributions rms radii. The fact that only 
q > 1.7 {m~^ cross sections are affected by RPA correlations suggests that this value of the 
momentum transfer constitutes the threshold above which the depletion of the inner regions 
of the target becomes sizable. However, since the density matrix used as an input to our 
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FIG. 5: Differential cross sections o-{9)/afiuth for protons incident on ^^O, ^''Ca and ^^Ca. Compar- 
ison between experimental data (symbols) and OMP predictions based on correlated (solid curves) 
and uncorrelated (dashed curves) descriptions of ground state. Cross sections offset factors and 
proton incident energies are indicated on the figure. Data are taken from Ref. for ^^O and 
40 Ca and from Q for^^Ca. 

microscopic OMP calculations conveys much more complex nuclear structure information 
that the radial density alone, disentangling the effects of the RPA correlations on nucleon 
scattering is a much more difficult task than the analysis of the q dependence of the radial 
density. Therefore, unlike the case of electron scattering, such an analysis can at best provide 
qualitative insight into the actual sensitivity of nucleon scattering to the presence of RPA 
correlations in the one-body density matrix of the target. 

2. Double counting 

Further tests of the sensitivity of our scattering predictions to changes in matter distri- 
butions have been performed by ignoring the Ap double counting correction terms (see ()17|) 
and (HH)). Elastic scattering calculation results performed with (solid curve) and without 
(dashed curve) these correction terms are shown in Figd for 201 MeV protons incident 
on ^o^Pb. First, Fig|7| shows that including or ignoring the Ap correction produces non- 
negligible changes in the calculated scattering cross section. Moreover, except for a local 
improvement a.t 6 = 54° over those using Ap ^ (solid curve), the agreement between data 
and the OMP calculation with Ap = is worse all over the range 6 > 34°. Setting Ap 
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FIG. 6: Proton elastic scattering from ^''^Pb: differential cross sections cr/aRuth as functions of the 
momentum transfer q. For more details, see caption of Fig. El 

to leads to increasing the rms radii from (r^/j)^^^ = 5.504 fm (Ap ^ 0, see Table |l]) to 
(r^^j)^/^ = 5.517 fm (Ap = 0), a value falling apart from the experimental result (r'^h)^^'^ ~ 
5.503(7) fm (see Table H)). The ^osp]-, neutron and proton radial shapes calculated assuming 
Ap = (dotted curves) and Ap 7^ (full curves) are shown in the insert of Fig|7| The above 
discussion shows that the Ap double counting correction to the RPA density matrix should 
not be ignored in scattering calculations. 



3. Skyrme Hartree-Fock model 

In recent years, Skyrme Hartree-Fock models have been considered to assess the neutron 
rms radius in 208p|j j2g|_ Furthermore, various Skyrme force parameterizations have been 
tested in NA g-folding model calculations to discern which one provides the best represen- 
tation of the neutron density. As a result, it turns out that SkM* seems appropriate when 
combining analyses of electron and nucleon scattering data, g-folding model calculations 
with HF/SkM* as input have again been performed and compared with calculations based 
on the present correlated ground state densities. The comparison made for (p,p) scattering 
off 208p]^ 201 MeV is shown in Fig|7| where the dotted and solid curves are for results 
from the HF/SkM* and HF+RPA/DIS based OMPs, respectively. The dotted and solid 
curve overlap each other over most of the angular range, except perhaps for angles above 
50". This is not surprising since both HF/SkM* and HF+RPA/DIS structure calculations 
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FIG. 7: Differential cross sections CF{9)/aRuth for 201 MeV protons incident on ^'^^Pb. Comparison 
between experimental data (symbols) and OMP predictions based on correlated (solid curves), 
correlated without double counting corrections (dashed curves), and Hartrec-Fock SkM* (dotted 
curve) descriptions of ground state. The insert shows comparison between proton and neutron 
radial densities for correlated (solid curves) and correlated without double counting corrections 
(dashed curves) descriptions of ground state. 

provide nearly identical radial matter distributions and neutron skins for ^°^Pb. However 
this similarity conceals more fundamental differences: whereas the SkM* interaction was 
designed to reproduce the measured charge radii of many stable nuclei within the HP frame- 
work only (its parameters take care of correlations present in nuclear ground states in an 
effective way at the mean field level), the DIS interaction is designed not to include such 
correlation effects in its parameterization, so that correlations can be explicitly taken care 
of, in a detailed way, at a level that goes beyond that of the mean field approximation. 

V. CONCLUSIONS 

We present a comprehensive analysis of ground state structure properties of doubly- 
closed-shell nuclei, together with the impacts they have on the interpretation of nucleon 
elastic scattering observables within the Melbourne g-folding model. Long range correla- 
tions are treated in the self-consistent RPA theory implemented with the DIS force, and the 
longstanding problem relevant to double counting is solved to calculate local and non-local 
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densities. The theoretical framework which in the past proved successful in the interpreta- 
tion of electron scattering measurements is shown to be equally successful in the analyses 
of nucleon elastic scattering up to 201 MeV. All the measured differential cross-sections, 
analyzing powers and spin-rotation functions are well described, without any adjusted pa- 
rameter. Turning off RPA correlations (or not implementing them properly, i.e. without 
considering double counting corrections) negatively affects the agreement between experi- 
mental data and calculations, an effect which gets more and more sizeable as incident energy 
and momentum transfer increase. It seems plausible that the differences observed between 
predictions are strongly tied to differences between correlated and uncorrelated matter den- 
sities at the surface and also towards nuclear center. Finally, since in the RPA theory, 
the correlated ground state happens to be the vacuum on which excited states are built as 
quasi-bosons excitations, a framework is at hand for extending our g-folding model analyses 
from elastic scattering to inelastic scattering from low to high excitation energy levels. Work 
along this line is in progress. 
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APPENDIX A: DEFINITION. 

The Hartree-Fock solutions in the spherical case take the form 

{x\{nlj),m,T) = Rl^iry [x^/^(<7) ® //^(r) . (Al) 

The operator (^^^^.^^^ creates a particle in this state and its hermitian conjugate defines 
the destruction operator a(nij),m,T- It is convenient to define destruction operators a(^nij),m,T 
through the relation 

Indeed, with this definition, both a'lnij) m t ^{nij),m,T transform under rotations like 
the component m of an irreducible tensor of rank j . Consequently creation operators of 
particle-hole pairs of definite angular momentum are readily constructed with the usual rules 
for coupling two tensors: 



K,h)i''^ = \^kr ^ Hh),r] ^ = I] C-^iM at),m„r Hh),mu,r ■ (A3) 



TUp, ruh 



The parity "tt" of the particle-hole pair that we indicate explicitly is defined by: n — {—Y^ 
As we did for the fermions, we define operators A 



^(,,,)(7r, J, M) = (-)^-^^(p_,)(7r, J, -M) , (A4) 

which annihilate particle-hole pairs of angular momentum J and projection M. As a con- 
sequence, by mixing and A, the Bogolyubov transformation defines operators and 
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© which respectively create and annihilate collective modes of definite angular momentum 
and parity. Finally, let us also recall that we consider only neutron and proton particle-hole 
pairs and consequently Tp = Th- 



APPENDIX B: THE QUASI-BOSON VACUUM. 

The expression of the quasi-boson vacuum takes a very simple form in the so-called 
canonical representation defined as follows 

(ph) 

{ph) 

The transformation D is an orthogonal transformation which mixes separately the creation 
and destruction operators of the original particle-hole pairs (it is orthogonal because our 
Bogolyubov transformation is real). It is defined by solving the eigen- values problem 

y- FyTT.jy^.jl D^^i pI'-^D""'-! , (B2) 



I {ph),(j>'h') 

In this representation, the vacuum reads 

io)=n n^h^«' ' 

■K,J \ a / 

with 

Z^\Y, th^^'^ ^ S+(7r, J, M)S„(7r, J, M) = ^ ^ th^^"^J [S+(7r, J) ® B^(n, J)] ° . 

(B3) 



2 / J a / J — av ' ' — / — uc\-- ■) - ■) — / 2 

■K,J,a M ■K,J,a 



The angle ^^''^ is related to the eigenvalues p^'-^ through the relation: 



V 1 + Pa 

This form shows clearly that the ^^'"^ 's provide a direct measure of the correlations which 
are induced by the RPA modes. 



[1] J. M. Cavedon et al, Phys. Rev. Lett. 49, 978 (1982). 
[2] E. N. M. Quint et al, Phys. Rev. Lett. 57, 186 (1986). 



18 



E. N. M. Quint et al., Phys. Rev. Lett. 58, 1088 (1987). 

I. Sick and P. de Witts Huberts, Comments Nucl. Part. Phys. 20, 177 (1991). 
V. R. Pandharipande and S. Pieper, Nucl. Phys. A507, 167 (1990). 
S. Fantoni and V. R. Pandharipande, Nucl. Phys. A427, 473 (1984). 
L. Lapikas, Nucl. Phys. A553, 297 (1993). 

Steven C. Pieper, R. B. Wiringa, and V. R. Pandharipande, Phys. Rev. C46, 1745 (1992). 
J. H. Heisenberg and B. Mihaila, Phys. Rev. C59, 1440 (1999). 
I. Bobeldijk et al, Phys. Rev. Lett. 73, 2684 (1994). 
L Bobeldijk et al, Phys. Lett. B353, 32 (1995). 

V. R. Pandharipande, C. N. Papanicolas, and J. Wambach, Phys. Rev. Lett. 53, 1133 (1984). 
V. R. Pandharipande, I. Sick, and P. K. deWitt Huberts, Rev. Mod. Phys. 69, 981 (1997). 
R. Wiringa, Nucl. Phys. A631, 70c (1998). 

Steven C. Pieper, V. Pandharipande, R. Wiringa, and J. Carlson, Phys. Rev. C64, 014001 
(2001). 

R. Wiringa, Steven C. Pieper, J. Carlson, and V. Pandharipande, Phys. Rev. C62, 014001 

(2000). 

B. Pudliner, V. Pandharipande, J. Carlson, Steven C. Pieper, and R. Wiringa, Phys. Rev. 
C56, 1720 (1997). 

S. Fantoni and V. R. Pandharipande, Phys. Rev. C37, 1697 (1988). 
V. Bernard and N. Van Giai, Nucl. Phys. A348, 75 (1980). 

C. Mahaux and R. Sartor, Nucl. Phys. A493, 157 (1989). 
C. Mahaux and H. Ngo, Nucl. Phys. A431, 486 (1984). 
C. Mahaux and R. Sartor, Nucl. Phys. A546, 65 (1992). 

C. Mahaux and R. Sartor, Adv. Nucl. Phys. 20, 1 (1991). 
J. Decharge and L. Sips, Nucl. Phys. A407, 1 (1987). 

J. Decharge, L. Sips, and D. Gogny, Phys. Lett. B98, 229 (1981). 
H. Lenske and J. Wambach, Phys. Lett. B249, 377 (1990). 

D. Gogny, in Lectures Notes in Physics, edited by H. Arenhorel and Dreschel (Springer, Berlin, 
1979), vol. 108, p. 88. 

Z. Y. Ma and J. Wambach, Phys. Lett. B256, 1 (1991). 

F. Petrovich, Nucl. Phys. A354, 499c (1981). 

W. Love, in The (p,n) Reaction and the Nucleon-Nucleon Force, Colorado, 1979 (Plenum 
Press, New York, 1980), p. 23. 

K. Amos, P. J. Dortmans, H. V. von Geramb, S. Karataglidis, and J. Raynal, Adv. Nucl. 
Phys. 25, 275 (2000), and references therein. 

H. V. von Geramb, in Lectures Notes in Physics, edited by H. V. v. Geramb (Springer- Verlag, 

BerHn Heidelberg New York, 1978), vol. 89, p. 104. 

P. Navratil and B. R. Barrett, Phys. Rev. C 57, 3119 (1998). 

P. Schwandt, in AIP Conference Proceedings, edited by H. O. Meyer (Springer- Verlag, Amer- 
ican Institute of Physics, New York, 1983), vol. 97, p. 89. 

J. F. Berger, M. Girod, and D. Gogny, Comput. Phys. Commun. 63, 365 (1990), and references 
therein. 

J. Blaizot and D. Gogny, Nucl. Phys. A284, 429 (1977). 

J. Bartel, P. Quentin, M. Brack, C. Guet, and H.-B. Hakansoon, Nucl. Phys. A386, 79 (1982). 
R. Machleidt, K. Holinde, and C. Elster, Phys. Rep. 149, 1 (1987). 
J. Raynal, computer code DWBA98, 1998, (NEA 1209/05). 



19 



[40] K. Amos, L. Canton, G. Pisent, J. P. Svenne, and D. van der KnijfF, Nucl. Phys. A728, 65 

(2003). 

[41] L. Canton, G. Pisent, J. P. Svenne, D. van der Knijff, K. Amos, and S. Karataglidis, nucl- 

th/0409050 (2004) submitted to publication to Phys. Rev. Lett. 
[42] F. Villars, in Fundamentals in Nuclear Theory, Collision Theory, Lectures presented at an 

internationnal course, Trieste, 1966, organized by ICTP, Trieste (IAEA, Vienna, 1967), p. 

269. 

[43] J. Decharge and D. Gogny, Phys. Rev. C 21, 1568 (1980). 

[44] J. Decharge, M. Girod, D. Gogny, and B. Grammaticos, Nucl. Phys. A358, 203c (1981). 
[45] P. Ring and P. Schuck, The Nuclear Many-body Problem (Springer- Verlag, New York, 1980). 
[46] B. Frois, J. B. BcUicard, J. M. Cavcdon, M. Huct, P. Leconte, P. Ludeau, A. Nakada, Phan 

Zuan Ho, and I. Sick, Phys. Rev. Lett 38, 152 (1977). 
[47] V. E. Starodunsky and N. M. Hintz, Phys. Rev. C49, 2118 (1994). 

[48] H. De Vries, C. W. De Jager, and C. De Vries, At. Data Nucl. Data Tables 36, 495 (1987). 
[49] B. C. Clark, L. J. Kerr, and S. Hama, Phys. Rev. C67, 054605 (2003). 

[50] A. Trzcinska et al., Phys. Rev. Lett. 87, 082501 (2001). 
[51] M. Csatlos et al, Nucl. Phys. A719, 304c (2003). 

[52] L. N. Blumberg, E. E. G. A. van der Wende, A. Zucker, and R. H. Bassel, Phys. Rev. 142, 
812 (1966). 

[53] A. Nadasen, P. Schwandt, P. P. Singh, W. W. Jacobs, A. D. Bacher, P. T. Debevec, M. D. 

Koritchuck, and J. T. Meek, Phys. Rev. C23, 1023 (1981). 
[54] W. T. H. Oers, Huang Haw, N. E. Davison, A. Igemarsson, B. Stagerstrom, and G. Tibell, 

Phys. Rev. CIO, 307 (1974). 
[55] V. Comparat, R. Frascaria, N. Marty, M. Morlet, and A. Wilhs, Nucl. Phys. A221, 403 (1974). 
[56] C. B. Fulmer, J. B. Ball, A. Scott, and M. L. Whiten, Phys. Rev. 181, 1565 (1969). 
[57] P. G. Ross and N. S. Wah, Phys. Rev. B140, 1237 (1965). 

[58] H. Sakagushi, M. Nakamura, K. Hatakama, F. O. A. Goto. T. Noro, H. Sakamoto, H. Ogawa, 

and S. Kobayashi, Phys. Rev. C26, 944 (1982). 
[59] P. Schwandt, H. O. Meyer, W. W. Jacobs, A. D. Bacher, S. E. Vigdor, M. D. Koritchuck, and 

T. R. Donohue, Phys. Rev. C26, 55 (1982). 
[60] D. A. Hutcheon et al., In Polarization Phenomena in Nuclear Physics, Proceedings of The 

Fifth International Symposium on Polarization Phenomena in Nuclear Physics, Santa Fe, 

AIP Conf. Proc, No. 69, (AIP, New York, 1981), p. 454 (edited by G.G. Holsen, R.H. Brown, 

N. Jarnie, W.W. Mc Naughton, and G.M. Hale, 1980). 
[61] M. Ju, Master's thesis, Simon Eraser University (1987). 
[62] N. Ottenstein, S. J. Wallace, and J. A. Tijon, Phys. Rev. C38, 2272 (1988). 
[63] L. Lee et al, Phys. Lett. B205, 219 (1988). 

[64] P. K. Deb, K. Amos, and S. Karataghdis, Phys. Rev. C 62, 037601 (2000). 

[65] J. H. Osborne et al, Phys. Rev. C70, 054613 (2004). 

[66] J. Klug et al., Phys. Rev. C67, 031601 (2003). 

[67] J. Klug et al., Phys. Rev. C68, 064605 (2003). 

[68] H. Seifert et al., Phys. Rev. C47, 1615 (1993). 

[69] A. E. Feldman et al, Phys. Rev. C49, 2068 (1994). 

[70] S. Karataghdis, K. Amos, B. A. Brown, and P. K. Deb, Phys. Rev. C 65, 044306 (2002). 



20 



